Checking 2 x M separability via semidefinite programming 
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In this paper we propose a sequence of tests which gives a definitive test for checking 2 x M 
separability. The test is definitive in the sense that each test corresponds to checking membership in a 
cone, and that the closure of the union of all these cones consists exactly of all 2xM separable states. 
Membership in each single cone may be checked via semidefinite programming, and is thus a tractable 
problem. This sequential test comes about by considering the dual problem, the characterization of 
all positive maps acting C 2x2 — * C . The latter in turn is solved by characterizing all positive 
quadratic matrix polynomials in a complex variable. 
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I. INTRODUCTION 

In the last decade entanglement has been recognized 
to be a fundamental notion in quantum information pro- 
cessing. While initially entanglement lead to "paradoxes" 
in quantum mechanics, it was later discovered to be a use- 
ful tool in teleportation, secure key distribution (quan- 
tum cryptography), quantum computation, etc. (see, 
e.g., and Q)- Despite its importance there does not 
yet exist a tractable conclusive test to check for entangle- 
ment or the lack thereof (separability). There are several 
partial tests the most famous of which is the Peres test Q 
which says that a separable state remains positive under 
taking partial transposes. In some cases the Peres test 
is also conclusive, namely in the low dimensional cases 
2 x M, M = 2, 3 (see Q) and in cases where the density 
matrix is of low rank (see Recently, in |(J the partial 
positive test was combined with an extendibility property 
to obtain a sequence of tests stronger than the Peres test, 
that all separable matrices need to pass. It is not known 
whether failure to pass this sequence of tests yields en- 
tanglement. Our main result concerns a sequence of tests 
that for the 2 x M separability problem do the opposite. 
Namely, passing one of the tests yields 2 x M separabil- 
ity This sequence of tests yields a definite answer to the 
2 x M separability problem in the sense that any 2 x M 
separable matrix is arbitrarily close to a matrix satisfying 
one of the tests. Let us present the details. 

Let M p (K) denote pxp matrices whose entries belong 
to K. Recall that a matrix p £ M2(Mu(C)) is called 



2 x M separable if there exists Xi 
that 



and yi e C M so 



Here t denotes the conjugate transpose and ® denotes 
the Kronecker product. Let n > 1 and p = I „^ ^ 



where A, B, and D are M x M matrices. We intro- 
duce the following convex set, which may be empty. 
Let Q{p;n) C M n+1 (M 2 {M M (€■))) consist of all posi- 
tive semidefinite block matrices T = (Tij)Zj=o (notation: 
r > 0) where 



r = 



,(a) 



eM 2 (M M (C)),0 <i,j <n, 



satisfy the conditions 
A, 



(i) E?=o r l a) 

(") E?= ri d) 



(iii) E"" 1 1 



■ii=0 2 



(T (c) 



= B, 



(iv) £?=o(r. 



(v) E£)*(r&*-r 



.(c) 
■ i,i 



r 4 ( 3) = 0, and 



i,i+kl 



= 0, k = 1, 



We will be interested in matrices p for which Q{p; n) is 
nonempty. Therefore we define 



A>, 



{ P eM 2 (M M (C)):g(p;n)^9}. 



It is straightforward to check that A n is a closed convex 
cone. Moreover, since Q{p;n) © 2 m Q G(p\n + \) one 
obtains that A n C A n +i, n> 1. 

It is an important feature of the set Q{p; n) that it is 
the intersection of the cone of positive semidefinite ma- 
trices (PSD) and an affine set. As such, the question 
of whether the set is empty or not, falls into a class of 
well studied problems, called semidefinite programming 
(SDP). Checking the nonemptyness of such an intersec- 
tion is called the feasibility problem in SDP, and sev- 
eral software packages (e.g., 0, 0) are readily avail- 
able to solve the feasibility problem numerically. A good 
overview article on SDP is 

We now state the main result. 



Theorem 1.1 The matrix p 
separable if and only if 



S M 2 (M M (C)) is 2 x M 
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p e \J n >iA n . 
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Based on the above theorem one may now formulate 
the following algorithm for checking for 2 x M separabil- 
ity. Given a time limit in which the check needs to be 
performed, make n as large as possible so that p G A n 
may be checked via SDP within the given time limit. In 
case the test comes out affirmatively, the given matrix p 
is 2 x M separable. When the test shows that p ^ A n the 
test is inconclusive. Still, given the content of Theorem 
11.11 the negative outcome may be an indication that the 
matrix is entangled. Of course, the larger n was chosen, 
the stronger the indication is. 

This is certainly not the first instance where SDP has 
been used for a quantum information problem. In fact, 
the earlier mentioned tests in Jfjl may be done by SDP. 
Moreover, in the papers 

[13 . E3» El> connections 
have been made between several other problems in quan- 
tum information and the versatile SDP tool. 

We will prove Theorem 11.11 by characterizing positive 
maps (also known as entanglement witnesses) acting on 
2x2 matrices. This characterization shall come about 
from the analysis of quadratic matrix polynomials in a 
complex variable. 



II. POSITIVE MAPS 



-<2x2 



iM KM 



We say that a linear map $ : — > n_/"~ j " is pos- 
itive when ^(G) > whenever G > 0. By linearity it 
suffices to check this condition for rank 1 matrices G. 
Furthermore, we may assume that G has the form 



G 



1 




or G 



If we let P = $(£u), Q = $(£a 2 ) and R = ${E 22 ), 
where Eij is the 2x2 matrix with a 1 in entry (i,j) and 
zeros elsewhere, checking $(G) > for all G as above 
amounts to checking that P > and 



\z\ 2 P + zQ + zQi + R > 0, z G 



(1) 



Since P > is automatically satisfied when (JJJ holds 
(let \z\ — ► oo ) it therefore suffices to check {Q. In other 
words, we have the following lemma. 

Lemma II. 1 The map $ : C 2x2 — > C MxM is positive if 
and only if P = $(J5u), Q = ${E 12 ) and R = ®(E 22 ) 
satisfy (JTJ- 

We therefore need to study quadratic matrix polyno- 
mials in a complex variable. 



III. QUADRATIC MATRIX POLYNOMIALS IN 
COMPLEX VARIABLES 



For M x M matrices P, Q and R, consider the matrix 
inequalities QJ. As observed, (JTJ implies P > 0, and 
clearly it also implies R > 0. Our analysis is based on the 



observation that we can eliminate argz and subsequently 
1 2 1 in inequality QJ. We do this as follows. Write z = 
re- 7 with r, 9 s R. Then 0) is equivalent to 

r 2 P + r(e^Q + e-^Q t ) +i? > 0, r,SeR. (2) 

We now use the following well-known fact. For complete- 
ness we provide a proof. Useful references on Toeplitz 
operators include 0, Chapter XXIII] and |l5j . 

Lemma III.l Given is the operator-valued trigonomet- 
ric polynomial H(z) := z^ 1 S^ + T + zS*, with S and 
T bounded linear operators on a Hilbert space Ti. Then 
H(z) > 0, \z\ = 1, if and only if the Toeplitz operator 



It s*t o o 

S T St o 
S T St 



(3) 



is positive semidefinite. 

Proof. Consider the multiplication operator g — > Hg 
on the Lebesgue space L 2 (T, 7i) of square integrable 
Lebesgue measurable functions on the unit circle T with 
values in 7i. By identifying L 2 (T,7{) with the Hilbert 
space £ 2 CZ,H) of square summable (in norm) sequences 
(^j')^=-oo) e the positive semidefiniteness of H(z) 
for all z G T is equivalent to the positive semidefinite- 
ness of the doubly infinite Toeplitz matrix A with symbol 
H(z) (that is, A is the doubly infinite version of ©)■ As 
P|) is the restriction of this doubly infinite Toeplitz ma- 
trix to £ 2 (N, H), positive semidefiniteness of Q follows. 

Conversely, let h — {hj)^^^ E £ 2 (Z,H) be so that 
hj = 0, j > B. Then it follows from the positive semidef- 
initeness of (0 that (Ah, h) > 0. Since sequences h of 
the above form are dense in £ 2 (Z, Ti), A > follows. But 
since H(z) is the symbol of this multiplication operator 
A, it follows that H{z) > 0, \z\ = 1. □ 

Fixing r and applying Lemma IIII.ll to ? we obtain 
that (0) holds for all 9 if and only if the following infinite 
block Toeplitz matrix is positive semidefinite: 



(R + r 2 P rQt -A 

rQ R + r 2 P rQt •• 
rQ R + r 2 P rQ^ ■■ 



\ 



> 0,r G 



7 



Notice that we may write (@J as 

r 2 A' + rL + N > 0, r G R, 



(4) 
(5) 



where K, L and N are infinite self-adjoint block Toeplitz 
matrices. Theorem 6.7 in 0] states that an oper- 
ator polynomial of even degree, 2q say, that is posi- 
tive semidefinite for all r G K allows a factorization 
T(r)T(ry , where T(r) is an operator polynomial of de- 
gree q. In our case we obtain that 

r 2 K + rL + N = (rT + S){rT + S) f . 
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This yields that K = TT\ L = TS^ + ST^ and N = SSl 
Consequently, if we let X = TS< — -j we get that X is 
skew-adjoint (X = —X^) and 



Applying this to our particular choice of K, L and N in 
terms of P, Q and R, we get that (Q holds if and only if 
there exists a skew-adjoint X — pQj)i,j>o so that 



K 



(6) 



Thus (JSJl implies the existence of X = —X^ such that © 
holds. The converse is also valid, as 



r-K + rL X=(rl >Ul K x 2 ^ X 



( 




v-t Q f . v-t 

"^00 2 "r^-lO 



Q i yt v-t Q i v-t 

2 ~r "^01 ^11 2 ' 21 



vv 



X, 



) 



9- + X±q Xn -5- + X12 • 



<R ••• 
R ■■■ 



"V 

/ 

/ 



> 0. 



(7) 



Performing a permutation, we may rewrite Q as 



A := 



P Aoo 
R 



X Q0 



xl 




Xqi 



p x u 

Xl R 




> 0. 



(8) 



Next we want to show that we may in fact choose Xij — 
Xi-j for all i and j. Perhaps the quickest way to do 
this is by using Banach limits (for the definition, see 0, 
Section III.7]). Carrying out some of the ideas related 
to Banach limits directly to the current situation, we 
obtain the following argument. Let A n be the infinite 
principal submatrix obtained from Ao by omitting the 
first n block rows and columns, each of size 2M. Thus 

A n has ^R^j m ^ e r ^S^ corner - Of course, 

A„ > for all n, and the sequence {A n } n >o is bounded 
(by ||Ao||). Consider now the bounded sequence {S„}„>o 
of averages defined via 



ent. Therefore, its limit Soo, which is positive semidef- 
inite, must have the same form as Ap as well, with Xij 



replaced by Yy, say. We claim that 



Y, 



for all i and j. Indeed, since limfc- 



lim*. 



n k + l 



= and 



n k +l 



= 0, we get that 



Y 



Y i+1 



lim — — 



X, 



k — >oc 



n k + 1 



lim 

k — >oc 



Xi+i j+i + • • • + X 



i+n k + l,j+n k +l 



Tl k + 1 



1 



(A H h A„) > 0, n G N . 



This sequence has a convergent subsequence {S nfc }fc>o 
in the weak operator topology. Notice that S nfc has 
the same form as Aq; only the operators Xy are differ- 



lim X i+n k +l,j+n k + l 



k^oo n k -\- 1 



n k + 1 



0. 



In conclusion, we have obtained that holds if and only 
if there exist Xj = —X_ it i G N so that 
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-Xq r 





Q Y 



x 1 







P X 

-X R 





-X_ 2 



*-2 









*_1 



> 0. 



(9) 



Let us denote the principal (n + 1) x (n + 1) block 
submatrix of @ by V(X , X n ; P, Q, R). Then we 
obtain the following result. 

Proposition III. 2 The matrices P, Q and R satisfy Q 
if and only if there exist Xi = —X\_ i , i = 0, 1, . . . so that 
V(X , ...,X n ;P,Q,R)>0 for all n. 



Let C n denote the cone of all matrices a = 



P Q 
Q f R 



so that there exist Xi = —X_ it i = 0,...,n, so that 
T>(Xq, . . . , X n ; P, Q, R) > 0. Using this notation we ob- 
tain the following. 

Proposition III. 3 The matrices P, Q and R satisfy Q 
if and only if 

a = ((J r) e n »>i C «- 



Combining Lemma III. II and Proposition IIII.3I we now 

obtain the following description of positive maps acting 

rn2x2 . rnMxM 



Proposition III.4 The map $ : C 2x2 -> C MxM is pos- 
itive if and only if 

\<P(E 21 ^(E 22 )J fcM »>i c «- 

Next we need duality to get back to the separability 
problem. We first prove the following auxiliary result. 

Lemma III. 5 The cones A n and C n are one another's 
dual. 

Proof. Let p G A n - In order to show that p G C* we 
need to show that trace(pcr) > for all cr e C„. Thus, 
let a G C n , and let Xi = —X_^ i = 0,...,n, be so 
that D := V(X , X n ; P, Q, R) > 0. Since p G A n 
there exists a V G Q(p;n). It is now straightforward to 
check that trace(pcr) = trace(rD), and since D and T are 
positive semidefinite, trace(per) > follows. This shows 
that A n <ZC*. 

For the converse inclusion, observe that when p = 

d) ^ ~^ n ' a ^ ^ sat ' s fy m S (i)"( v ) we 



have that (rjj)™-_ 0. Notice that the block matri- 
ces (ry)™-_ satisfying (i)-(v) describe a afiine space, 
which we may denote as G + C where G is a fixed ma- 
trix and £ is a linear subspace (described by all ma- 
trices satisfying (i)-(v) with A = D = B = 0). This 
affine space G + C is separated from the cone PSD of 
positive semidefinites. Thus by the Hahn-Banach the- 
orem and the selfduality of PSD there exist a positive 
semidefinite W so that trace(M / (G + L)) < for all 
L G C. But then W G C 1 - and thus W is of the form 



V(X ,...,X n ;P,Q,R). Let a = 



P Q 

Q f R 



which be- 



longs to C n . As trace(cr j o) = trace(W / G) < 0, it follows 

that p g C*. Thus A n = C*. 

Since C„ is closed, A* n = C** = C n follows also. □ 
We are now ready to prove the main result. 

Proof of Theorem II. II By |4| we have that the cone of 

2 x M separable matrices has the dual 



{ 



f$(E n ) $(£ 12 ; 
\${e 21 $(£ 22 ; 



^2x2 



is positive }, 



which by Proposition 1111.41 equals H^ =1 C n - The cone of 
2 x M separable matrices therefore equals 



In the last step we used Lemma fill. 51 □ 
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